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Abstract. We devise a scalable scheme for simulating a quantum phase transition 
from paramagnetism to frustrated magnetism in a superconducting flux-qubit network, 
and we show how to characterize this system experimentally both macroscopically and 
microscopically. The proposed macroscopic characterization of the quantum phase 
transition is based on the transition of the probability distribution for the spin-network 
net magnetic moment with this transition quantified by the difference between the 
Kullback-Leibler divergences of the distributions corresponding to the paramagnetic 
and frustrated magnetic phases with respect to the probability distribution at a 
given time during the transition. Microscopic characterization of the quantum 
phase transition is performed using the standard local-entanglement-witness approach. 
Simultaneous macro and micro characterizations of quantum phase transitions would 
serve to verify a quantum phase transition in two ways especially in the quantum realm 
for the classically intractable case of frustrated quantum magnetism. 
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1. Introduction 

Experimental quantum simulation [1] opens vistas for exploring foundational quantum 
principles such as studies of Quantum Phase Transitions (QPT) [2-6]. Although 
the phase of matter can only be defined macroscopically, the QPTs are currently 
explored through microscopic characterizations that depend on local observables such as 
individual-particle spin states or few-body states or quantum entanglement witnesses [7] . 
Macroscopic probing, which measures some global property of a system without having 
any access to individual particles, is not currently employed in quantum-phase-transition 
experiments. Yet conducting joint macro and micro characterizations of quantum phases 
and transitions between them is important for self-consistent verification of QPTs in the 
macro and micro regimes, assuming that one may not necessarily imply another [8,9]. 

Here we propose a scalable scheme comprising nearest- and next-nearest-neighbor- 
coupled spin-half network of radio-frequency superconducting-quantum-interference- 
device (rf-SQUID) flux qubits [10, 11], where the two levels comprise one clockwise 
and one counter-clockwise super-currents. The clockwise super-current state |0) is 
equivalent to a spin-down flux state IJ,), and the counter-clockwise super-current 
state |0) corresponds to a spin-up flux state |t)- This flux-qubit spin network will realize, 
and admit macro and micro characterization of, a controllable QPT from paramagnetism 
to frustrated magnetism [12]. 

We focus on quantum simulation of the paramagnetic to frustrated-magnetic 
QPT [6] because of the need to simulate the fascinating behavior manifested in frustrated 
spin systems [12], such as spin glasses, spin ices and quantum spin liquids due to 
competing spin-spin interactions. These competitions give rise to an exponentially large 
subspace of degenerate ground states for such frustrated networks [13], which strongly 
motivates quantum simulation of these systems [1]. We seek to identify the macro 
signature of the onset of quantum frustrated magnetism from paramagnetism through 
the statistics of the order parameter rather than solely by its mean. 

Quantum simulations of frustrated Ising networks [14] have been performed with 
ion traps [7], quantum dots [15], and nuclear magnetic resonance [16] and has been 
proposed for Rydberg atoms [17]. Larger frustrated systems have been realized for 
optical systems [18] and trapped ions [19]. All these schemes employ microscopic probes 
of frustration and are not yet amenable to macroscopic probing, which is why we propose 
a new setup corresponding to a scalable flux-qubit superconducting architecture that 
admits probing of both macro and micro signatures of frustration. 

The rest of the paper is organized as follows: Section 2 describes the physical model 
for the architecture considered in this work. In Section 3 we elaborate the control 
procedure and discuss how micro and macro signatures can be extracted for our system. 
Section 4 shows the results of our simulation and we conclude in Section 5. 
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Figure 1. Schematic for (a) n nearest-neighbor- and next-nearest-neighbor-coupled 
qubits and (b) three coupled qubits in a triangular architecture, which is the smallest 
lattice required for frustrated magnetism. Rounded squares denote qubits, and solid 
lines denote couplings. 


2. Physical model 

Our goal is to analyze the model for nearest- and next-nearest-neighbor-coupled rf- 
SQUID superconducting flux qubits (as shown in Fig. la) and devise a scheme to 
extract the macroscopic as well as microscopic signatures of QPT from paramagnetism 
to frustrated magnetism. However, we hrst consider three rf-SQUID superconducting 
flux qubits coupled to each other via tunable inductive couplings (shown in Fig. lb), 
as this triangular architecture is the smallest lattice that can achieve the regime of 
frustrated magnetism. 

2.1. CCJJ flux qubit 

The Compound-Compound Josephson-Junction (CCJJ) |0) shown in Fig. 2 is proposed 
here as our favored flux qubit (with basis states |0) and |c!))) due to its exquisite 
tunability [10,11] and its demonstrable scalability in the sense that hundreds of such flux 
qubits can be coupled in a spin network [20-22]. The two flux-qubit levels correspond 
to the clockwise and anti-clockwise circulating currents in the qubit loop (denoted ‘Q’ 
in Fig. 2). 
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Figure 2. A circuit diagram of an rf-SQUID CCJJ superconducting flux qubit, where 
the quantized circulating current Iq across the qubit loop (denoted by ‘Q’) defines the 
|0) and |0) states of the qubit. The CCJJ loop comprising four Josephson junctions 
is coupled to a bias current Jccjj for controlling the height of the flux-qubit potential- 
energy barrier. The time-dependent bias current /ccjj(f) induces a time-dependent 
magnetic flux <i)ccjj(f) across the CCJJ loop. This flux controls the |0) O |0) 
tunneling probability. The compensator C, which is biased by current Iq via an in 
situ tunable inductive coupler to the flux qubit, is controlled by the static flux bias 
$c- It compensates for the changing current Iq, such that the qubit frequency e can be 
maintained to a specific value during the experiment. The readout circuitry (denoted 
‘RO’) can detect any variation of the circulating current Iq and is used to perform the 
readout operation on the flux qubit. 


In the {|0), lO)} basis, the Hamiltonian of the CCJJ flux qubit is off-diagonal due 
to the tunneling between these states across the hnite potential barrier. We, therefore, 
have two independent control parameters: e that controls the asymmetry between two 
potential wells, and A that controls the potential barrier height between two wells. The 
Hamiltonian of this qubit is [10,11] Hq = — {ea^ + Aa^) for e = |Jq| (<hQ — <ho) with Iq 
the super-current circulating in the rf-SQUID loop, <hQ the time-dependent qubit-flux 
bias applied across the qubit, and the three Pauli matrices. It is possible to turn 
off the diagonal terms of the Hamiltonian by setting <I)q = <I)o, a regime often referred 
to as the degeneracy point. Tunneling energy A is a function of ‘hccjj that can be 
controlled via Jccjj- 

Fabricating such qubits having identical controllability but superior energy 
relaxation time (expected Ti > 10 /rs) is underway [23], which could serve as an 
ideal platform for conducting such QPT experiments. Coherence times of tens of 
microseconds are orders of magnitude higher than our operation time for simulating 
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the QPT, which guarantees that the effect of decoherence-induced noise is negligible for 
our case. On the other hand, superconducting qubits operate at around 10 mK [23-25], 
which corresponds to a thermal excitation with frequency ~ 1.3 GHz, which is far 
detuned from the frequencies of the quasiparticles in the flux qubit (~ 6 GHz.) [21]. 
This detuning implies that the mean density of thermal photons in the system at the 
operating temperature is less than 1%, so the possibility of absorbing thermal excitations 
by the system is minimal, thereby rendering our architecture resilient against possible 
environment-induced noises. 



Figure 3. A circuit diagram of three rf-SQUID flux qubits coupled with each other 
via inductive tunable couplers. The square around the system couples all three flux 
qubits and measure the total magnetic moment of the system which we use towards 
the macroscopic characterization of the QPT in the network. 


2.2. Triangular architecture of coupled CCJJ qubits 

We now consider a triple-spin network with the three GGJJ flux qubits inductively 
coupled in a triangular arrangement shown in Fig. 3. In this architecture, each pair of 
qubits is coupled via an inductive coupler that can be tuned in situ via the external 
flux biases [11]. Such an inductive coupling scheme induces a a-af-type interaction 
between the i*'^ and the qubits with an adjustable coupling strength Jij [11]. The 
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Hamiltonian for this network of three coupled flux qubits is 

3 2 3 

". = E (1) 

2=1 2=1 J =2 + l 

with i,j denoting qubit indices and Hq the Hamiltonian for the flux qubit. Using 
the control mechanisms described earlier, we can achieve e* = 0 by setting each external 
qubit flux bias at the degeneracy point of the qubit. 

3. Simulating QPT from paramagnetism to frustrated magnetism 

In this section we discuss how to control the quantum Hamiltonian (1) in order to achieve 
the frustrated regime starting from the paramagnetic phase. We compute the fidelity 
susceptibility for the ground state of the system and explore its divergent behavior, 
which in fact denotes a QPT without any prior knowledge of any local order parameter. 
Finally we discuss how micro and macro signatures of the QPT can be extracted for our 
system. 


3.1. Control scheme 

In our approach to simulating frustrated antiferromagnetism, we initially prepare the 
ground state of the system Hamiltonian (1) with e* = Jij ~ 0 Vi,j and then vary A 
and J (assuming Aj = A and Jij = J for all i,j) adiabatically via d’ccjj and 
respectively. The control pulses for these parameters are designed such that, at t = 0, 
J/A -C 1 (~ 10“®), and, at t = tfinah J/^ 3> 1 (~ 10®) [21,26]. The Hamiltonian 
ground state \g) at t = 0 is \g(t = 0)) = |+ + +) [27], which is separable. During the 
system’s adiabatic evolution towards the frustrated regime, the three qubits become 
entangled, and the system ground state at t = tfinai is (neglecting global phase) 

(un> + im> + itu) - im> - int> - im»/V6. 


3.2. Fidelity susceptibility 


Fidelity susceptibility offers a measure for QPT in absence of any prior knowledge of 
order parameters [28,29]. The fidelity susceptibility xf quantifies the drastic change in 
the ground state of the quantum system during the phase transition, and is defined as 


Xf(A) — 


a2ln|(To(A)|To(A + 6))| 




( 2 ) 


e=0 


where Tq is (normalized) ground state wavefunction of the system and A is the control 
parameter, which is J/A for our case. It is possible to show that Xf(A) > 0 for all A, 
and can also be expressed as [29] 



d^o\ 
dX / 



aTo\ /9To 
dX /\d\ 



( 3 ) 


We use Eq. (3) to compute the xf for our system, and seek its divergent behavior that 
denotes the QPT. 
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Figure 4. Plots showing the fidelity susceptibility yp as a function of time for n 
nearest- and next-nearest-neighbor-coupled flux qubits. The divergence in yp denotes 
a QPT. 


Fig. 4 shows the hdelity susceptibility xf as a function of time for various lattice 
sizes. During the time f = 0 to f = 50 ns, we simultaneously vary J from 0 to 5 GHz 
and A from 5 to 0 GHz adiabatically. We observe that while yp ^ 1 for most of the 
time during the change of our control parameter, it however shows divergence when 
J//S.^ 1. The plots in Fig. 4 not only indicates the existence of a QPT for our system, 
but also shows that extracting the signatures of QPT should be possible for a lattice 
with n ~ 12 for which the divergence of yp is prominent. This requirement is especially 
useful as the classical simulation of any quantum spin system becomes exponentially 
expensive with higher values of n. For the macroscopic characterization of QPT, we 
therefore restrict ourselves to a lattice with n = 12, which turns out to be sufficient for 
the architecture considered in this work. 

3.3. Microscopic signature 

QPTs are usually characterized microscopically by the divergence of correlation length 
near the critical point [2]. Experimentally, however, it is not always efficient to determine 
this divergence, and, therefore, alternative system-specihc schemes are often used, such 
as measurement of some local entanglement-witness operator for spin-systems [7]. 

The QPT for our architecture can be characterized by selecting a block of 
three adjacent flux qubits (which is in fact a triangular architecture as shown in 
Fig. lb) and measuring a three-qubit entanglement-witness operator on that block. 
Such a measurement is performed locally on a few-spin subsystem of the entire 
lattice. Moreover, the measurement of entanglement-wit ness operator requires readout 
operations on each individual qubit, hence is a microscopic characterization. 

The symmetric W-state witness operator [7,30] 


















satisfies (H^s) > 0 for separable states and (ffs) < 0 for entangled states. For onr case, 
we expect (fFs) = — 2 at t = 0 and (PFs) = — 3 at t = ffinai, which prominently 

signihes the entanglement change. 

3 . 4 . Macroscopic signature 

For macroscopic characterization, we can place a dc-SQUID loop aronnd the entire 
system, snch that it couples all the flux qubits (Fig. 3 shows such a loop for the three- 
qubit case) and measure the total magnetic moment fi along the ^-axis [31]. The net 
magnetic moment of the entire system along z-axis can be expressed as a sum of the 
z-components of the individual spin magnetic moments entire system. 

Classically, the total magnetic moment of a system per unit volume is the magnetization 
of the system, which is a macroscopic quantity. We refer to this readout procedure as a 
macroscopic characterization of frustration, as it does not measure any local properties 
of the system addressing each flux qubit individually (elaborated in Appendix A). 

The spectrum {/i^} of the total magnetic moment operator assumes only integer 
values. Repeated measurements yield a probability distribution for {/i^} centered 
at = 0 for both the frustrated as well as paramagnetic phases. The tail of the 
distribution differs between frustrated and paramagnetic phases. For the paramagnetic 
phase, the state of the entire system is which indicates that all states are equally 

probable, if measured in the {|t), li)}®"' basis. Therefore, the probability that I spins 
are in |)") state and I -|- k spins (with n = 21 + k) are in |J,) state (for which \p,^\ = fc), 
^((n-fc)/ 2 ) ((n+fc)/ 2 )’ ^hich meaus that the probability distribution of /i^ is a binomial 
function for the paramagnetic case. For the frustrated phase, the probability distribution 
depends on the geometry of the lattice, which is why it is called geometric frustration. 
Whereas, an exact closed-form analytic solution of the distribution for this case is hard 
to derive, the frustrated antiferromagnetism implies that the peak of the distribution is 
still centered at = 0, and the tail is an exponential decay (as opposed to binomial), 
which can be verihed by numerical simulations. 

In order to characterize the change in the tail of the probability distribution, we 
introduce a measure 

V[V{t)] := Z^KL[Pexp,P(t)] - /^KL[Pbin,P(t)], (5) 

on the probability distribution 'P{t) at time t with Pexp and Pbin best-£t exponential 
and binomial distributions corresponding to the hnal and initial states. The Kullback- 
Leibler divergence Pkl := J2k^k^og{Pk/Qk) quantihes distinguishability between two 
probability distributions P and Q. If the distribution V{t) changes from binomial to 
exponential, then P goes from positive to negative, which indicates a phase transition 
from paramagnetic to frustrated magnetic phase has occurred. 
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Figure 5. Plots showing the microscopic signatures of QPT. (a) The control pulse 
where A gets turned off and J turned on adiabatically with time, (b) The probability V 
for each state (denoted by subscripts in the legend) with time under the control pulse, 
(c) Expectation of the symmetric IT-state witness operator as a function of time under 
the control pulse. The change of sign indicates the QPT. 


4. Results 

Figure 5 shows the results for the microscopic characterization of the QPT for our three- 
qubit device. The adiabatic pulses for A and J are shown in Fig. 5(a) for couplings 
slowly turned on and local controls off. Figure 5(b) shows the probability distribution for 
all eight states. Whereas at t = 0 the probabilities are all equal, they gradually decay 
for Ittt) and im) states with a uniform probability distribution for the remaining 
states, which is expected for a paramagnetic to frustrated magnetic phase transition. 
Figure 5(c) shows the change of symmetric hF-state entanglement-witness operator with 
time. At t = tfinai the state of the system is a superposition of two W states, and therefore 
(IF 5 ) goes from positive to negative indicating the generation of entanglement in the 
system, which we consider as our microscopic signature. 

Figure 6 shows the macroscopic signatures of the QPT for an n-qubit network of 
nearest-neighbor- and next-nearest-neighbor-coupled superconducting flux qubits shown 
in Fig. la. We plot the initial and hnal probability distributions for the total magnetic 
moment of the system in Fig. 6 a for n = 12. Although for both cases the distribution 
is centered around the origin, a signihcant difference is observed in the tails. The 
tail of the probability distribution changes from a binomial to an exponential decay. 
'D\P{t)\ quantihes this change and 6 b shows how V changes gradually from positive to a 
negative value with time signifying a QPT with respect to a macroscopic probe, which 
is the SQUID loop around the network for our case. 
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(a) 



(b) 



Figure 6. Macroscopic signatures of the QPT. (a) The probability distributions V{p,^) 
of the total magnetic moment for the network of 12 flux qubits at initial and final 
times. The black dots denote the actual data computed via numerical simulation and 
solid lines denote interpolations through the data. The tail changes from binomial 
to exponential denoting the QPT. (b) Our distance measure based on the Kullback- 
Leibler divergence as defined in Eq. (5) as a function of time. The change in sign 
indicates the QPT. 


5. Conclusions 

In this work, we have put forward a realistic scheme to simulate a QPT, where a network 
of nearest-neighbor- and next-nearest-neighbor-coupled superconducting flux qubits is 
subjected to a transition from the paramagnetic phase to the frustrated magnetic phase. 
Characterizing a QPT with a micro as well as a macro probe is a conceptual requirement, 
rather than an experimental preference. 

Whereas the microscopic characterization of the phases can be performed measuring 
the appropriate entanglement witness, we have introduced a measure in this work based 
on the Kullback-Leibler divergence that can characterize the QPT macroscopically. 
Our proposal for such simultaneous characterization of QPTs at macro and micro 
scales enables a rigorous authentication of various quantum phases, as well as offers 
a standardized procedure to self-consistently verify the signatures of quantum phase 
transitions in upcoming experiments. 
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Appendix A. Kullback—Leibler Divergence: A Macroscopic Measure for 
QPT 

A phase transition, classical or quantum, is always dehned for a macroscopic system. 
Here we are interested in the question if the phase transition can be characterized 
macroscopically or not; in other words, if the probe for the QPT is macroscopic or 
not. The motivation for this question is described in the previous section, and here 
we elaborate on why the SQUID loop considered for characterizing the transition from 
paramagnetic to frustrated magnetic phase is in fact macroscopic. 

The SQUID loop around the entire network of nearest- and next-nearest-neighbor- 
coupled flux qubits is measuring the total magnetic moment of the system, which in 
terms of the individual spins can be expressed as, However the SQUID 

loop does not resort to measuring each individual spins and then summing it up; in fact 
it does not even have any access to each individual spins in the network. Therefore, this 
probe functions as a macroscopic probe. 

Since the spins in the network have many possible arrangements, the SQUID loop 
can End different values for the total magnetic moment of the system corresponding 
to various spin-arrangements, and by repetitive measurements one can construct the 
probability distribution for these different arrangements. As the control parameters are 
varied from the paramagnetic regime to the frustrated magnetic regime, the probability 
distribution also changes, and the question of characterizing the phase transition boils 
down to the question of how to distinguish between two probability distributions 
quantitatively. 

The probability distributions corresponding to paramagnetic and frustrated 
magnetic phases are shown in Fig. 5. The peaks for both of these probability 
distributions are centered about the origin (|/i^| = 0), whereas their tails differ 
signihcantly. We thus point out that the tail of the probability distribution carries 
the information about the phase. In order to quantify the change in the tail of 
the distribution we propose to compute the Kullback-Leibler divergence Dkl, a non- 
symmetric measure of the distance between two probabilistic distributions. It is 
important to mention in this context that the measures like Kullback-Leibler divergence 
or Jensen-Shannon divergence has already been used to characterize the classical phase 
transition (not as a macroscopic measure) from Bose gas to Bose-Einstein condensate, 
where the transition is also characterized by the change in the tail of the velocity- 
distribution [32]. 
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The measure we introduced here to characterize the QPT (given by Eq. (3)) takes 
positive value for paramagnetic and negative value for frustrated magnetic phase, and 
thereby quantihes not only the change in phase, but also indicates the phase itself 
for the entire system. It is important to note that the Kullback-Leibler divergence 
is a measure of classical relative entropy, whereas characterizing a QPT for our case. 
This is not surprising, since the change in the probability distribution for different 
phases are still mediated by quantum fluctuations, and a classical measure is sufficient 
to quantify this change. One can possibly introduce some macroscopic measure of 
entanglement contained in the entire network, and characterize the QPT based on that 
measure. However, we do not attempt to dehne such a measure as it does not change 
the conclusions of our work. 
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